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ABSTRACT 

This is the second article of a series devoted to European Extremely Large Tele- 
scope (E-ELT) site characterization. In this article we present the main properties of 
the parameters involved in high angular resolution observations from the data collected 
in the site testing campaign of the E-ELT during the Design Study (DS) phase. Ob- 
servations were made in 2008 and 2009, in the four sites selected to shelter the future 
E-ELT (characterized under the ELT-DS contract): Aklim mountain in Morocco, Ob- 
servatorio del Roque de los Muchachos (ORM) in Spain, Macon range in Argentina, 
and Cerro Ventarrones in Chile. The same techniques, instruments and acquisition 
procedures were taken on each site. A Multiple Aperture Scintillation Sensor (MASS) 
and a Differential Image Motion Monitor (DIMM) were installed at each site. Global 
statistics of the integrated seeing, the free atmosphere seeing, the boundary layer see- 
ing and the isoplanatic angle were studied for each site, and the results are presented 
here. In order to estimate other important parameters such as the coherence time of 
the wavefront and the overall parameter "coherence etendue" additional information 
of vertical profiles of the wind speed was needed. Data were retrieved from the Na- 
tional Oceanic and Atmospheric Administration (NOAA) archive. Ground wind speed 
was measured by Automatic Weather Stations (AWS). More aspects of the turbulence 
parameters such as their seasonal trend, their nightly evolution and their temporal sta- 
bility were also obtained and analyzed. 

Subject headings: site testing - site characterization 
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1. Introduction 



The site selection for the future large European telescope is a fundamental issue and was 
undertaken within the E-ELT Design Study proposal funded by the European Community. First 
meetings and contacts to define the site selection project started in 2003. Possible interested 
partners and institutions were approached and a first version of design and plans was submitted to 
the European Commission in February 2004. After revision, using the committee feedback, the 
final proposal was acce pted at the end of 20 04. The Site Selection work started formally in 2005 



Vemin et al. 



(|201ll) (hereafter Paper I) we presented an overview of the 



to end in May 2009. In 
campaign. 

In the present work the statistics of the parameters relevant to high angular resolution 
(HAR) astronomy of such a large telescope are presented and discussed in detail. Results were 
obtained after the analysis of about one y ear of atmospheric tur bulence observations with the 



same instrumentation (MASS-DIMM, see 



Kornilov et al, 



(l2007h ) in four different sites: Aklim (in 



Morocco), Macon (in Argentina), Observatorio del Roque de los Muchachos (ORM) (in Spain) 
and Ventarrones (in Chile). This study is similar to the site characterizat ion produced by the Thirty 



Meter Telescope (TMT) team: its first article fr omlSchock et a 



statistics of the seeing a nd the isoplanatic angle 



Skidmore et al. 



Travouillon et al, 



(120091) 



(20091) and those focused on the 



( 20091) . and on the coherence time 



The paper is organized as follows: First, the overall observing configuration at each site is 
detailed in Sec. II. H and the definitions of the parameters under study are introduced in Sec. 11.21 
The employed instruments and the data provided by them are described in Sec. 12.11 while Sec. 12.21 
is devoted to the archives used in order to obtain the complete vertical wind profiles at each site, 
which are needed to compute some of the parameters. Their main global statistics is covered in 
Sec. 121 The seasonal behavior, the evolution during an averaged observing night and the stability 
of the studied atmospheric parameters are addressed, respectively, in Sees. |4l[5]and[6l Finally, the 
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main results are summarized in Sec. U\ 



1.1. Observing configuration 



Our aim was to monitor the atmospheric turbulence at the four candidate sites using well- 
known, reliable and as much homogeneous instrumentation as possible. All the MASS -DIMM 
instruments were installed on 5 m high towers. This was ag reed based on previous studies on the 



surface layer thickness (see 



Vemin & Munoz-Tun6nl(ll992|) ). Here follows the instrument setup: 



o Four identical MASS-DIMM instruments. 



o Four identical telescopes Celestron CI 1 (11 inches). 



o Four identical fast read-out CCDs for DIMM devices: the PCO PixelFly VG 



o Each MASS device contains four photo-multipliers (|Kornilov et al.ll2OO70 



o Four towers in order to observe at 5 m above ground level. 

o Three robotic mounts (ASTELCO NTM-50(|^) installed at ORM, Macon and Ventarrones, 
and one automatic (Losmandy Geminj^ installed at Aklim. 

o Four Automatic Weather Stations (AWS) at few meters above ground level (see Paper I for 
details). 



The observations taken at Aklim and ORM sites were carried out regularly by observers (in 
particular, at ORM site observations took place five nights per week; at Aklim, observations are 



^http : //www ■ pco ■ de/sensitive- cameras/pixelf ly- vga/ 



jhttp : //www ■ astelco . com/product s/ntm/ntm . htm 



- jhttp : //www . losmandy . com/losmandygoto/gotospecs . html 
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less regular due to strong difficulties to access the mountainous site), while those taken at Macon 
and Ventarrones sites are obtained robotically every night. More details about the duty cycle of 
the MASS-DIMM instrument in each site was already given in Fig 1 of Paper 1. 

Concerning the observation itself, the configuration adopted in each site is summarized in 
Tables [Hand [21 The measurements taken with both instruments, DIMM and MASS, are filtered 
according to well defined criteria stated in Sec. 12.1.41 After filtering, one obtains the total number 
of accepted data N^^cc and the percentage of rejected data N^ej- These parameters, together with the 
total number of exposures per measurement, A^exp, the exposure time, fexp and the median sampling 
time. At are given in Tables [Hand [2] for the DIMM and the MASS respectively. 



The major parameters relevant to high angular resolution (imaging, adaptive optics, 
interferometry) have been grouped into two classes: "integrated" parameters and "profiles". The 
later class is represented by optical turbulence profiles, Cj^ih) and wind speed profiles V(h). It 
is well known that integrated parameters, such as seeing or Fried's radius, isoplanatic angle and 
coherence time, can be deduced from both the above-mentioned profiles as: 

• Fried's radius: 



1.2. Parameters 




(1) 



• Seeing: 




(2) 



• Isoplanatic angle: 




(3) 
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• Coherence time: 



To = O.OSSA^'^ 



( CO 

vo 



-3/5 



(4) 



where the light wavelength is X = 0.5 /um and all the measurements are referred to zero zenith 
angle. 

Although seeing is generally considered the most important parameter for HAR astronomy, 
various combinations of s, 6, r are also used to compute some figure of merit, as already discussed 
by Paper I, dependin g upon the hi gh angular resolution technique employed. A more general 



approach is given by 



Lloydl (|2004r) . who defines the "coherence etendue" Go, in which a photon 



remains coherent. Gq takes into account a combination of Fried's radius, isoplanatic angle and 
coherence time: 

Go = rlToel. (5) 

This new formulation shows a strong dependence of Go with ro and and less with To. Gq is 
computed with ro, tq and Oq respectively expressed in m^, ms and arcsec^. 



2. MASS-DIMM and complementary NOAA data 
2.1. MASS-DIMM 

MASS and DIMM devices are attached to the same equatorial mount and track at the same 
star, but each instrument has its own setup. 



2.1.1. Differential Image Motion Monitor (DIMM) 



DIMM provides accurate, absolute and reproducible integrated seeing data a 
systematic control tests on the focus or saturation are however important (see e.g. 



though 



Tokovinin 
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2002 



Varela et al 



2004). The instrument description is given in 



Sarazin & Roddier 



1990 and 



Vemin & Munoz-Tunon 



1995 



Since the early nineties, DIMMs have become very popular and 
have been used at different observatories. DIMMs are now auxiliary instruments for telescope 
operation and complement Adaptive Optics (AO) experiments. For what concerns the site 
selection, accurate statistics is an important issue. A lot of results have been recorded i n large 



databases that can be available in iMunoz-Tunon et al. 
example, in 



Ehgamberdiev et al, 



1997 : 



Ehgamberdiev et al. 



2000. For 



20001 seeing values at Paranal and La Palma were analyzed and 



compared for more than two years. From this analysis, the excellent behavior of the two sites is 
clear and reinforces their pre-selection for hosting the future E-ELT. But the DIMM provides only 
the seeing, and a C^(/z) profile is needed to access the isoplanatic angle (Eq. [3]). 



2.1.2. Multi-Aperture Scintillation Sensor (MASS) 



This instrument detects fast intensity variations of light in four concentric apertures using 
photo-multipliers. Every minute, the accumulated photon counts obtained with micro-exposures 
of 1 ms are converted to four normal scintillation indices and to six differential indices for each 
pair of apertures. This set of ten numbers is fitted by a model of six thin turbulent layers at 



pre-defined altitudes, hj = 0.5, 1, 2, 4, 8, and 16 km above the site altitude (|Kornilov et al 



2003 ) 



Another model of three layers at floating altitudes is fitted as well. The set of integrals of the 
refractive index structure constant. 



Ji 



/ 

Jm 



Cjf{h)dh, 



(6) 



(th layer 



in these six (or 



Tokovinin et al 



t hree) 



layers represent the optical turbulence profiles measured by MASS (see 



(120031) for details on MASS weighting functions). Turbulence near the ground 
does not produce any scintillation: MASS is blind to it and can only measure the seeing in the free 
atmosphere. 
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MASS has been cross-compared wit h the Generalized Scintillation Detection And Ranging 



(G-SCIDAR ) optical turbulence pr ofiler (|Avila et al 



Mauna Kea (Tokovinin et al, 



1997b . during a campaign performed at 



2005|) showing very good agreement. SCIDAR has proved to be 



the most efficient and reliable technique to accurately measure the optical vertical structure of 
the atmospheric turbulence strength from ground level, although it requires a one-meter class 
telescope to perforrn the observations. A more recent study carried out at Paranal Observatory 



(IDali Ali et al 



2010|) also produced consistent results. Similar comparisons between the 
parameters provided by the MASS-DIMM instrument and the G-SCIDAR, made at ORM, will 
be addressed in the forthcoming issue within the present series which will be devoted to the 
G-SCIDAR profiles. 

Assuming that the optical turbulence profile remains constant within each slab defined by the 
MASS, one can deduce the isoplanatic angle (Eq. The coherence time (Eq. H)) still requires the 
knowledge of the wind speed profile which is not delivered by the MASS, but will be retrieved 
from meteorological archives, as explained in Sec. 



2.1.3. Cross-calibration of the DIMM device 

This section describes the comparison made between the seeing values obtained with one 
of the DIMM devices employed during the E-ELT site characterization project (the DIMM part 
of the MASS-DIM M instrument) and with an exis ting stable seeing monitor at ORM (hereafter 



called lAC-DIMM, 



Vernin & Munoz-Tunon 



19951). An on-line report is availablqj with detailed 



information on the systems setup and on data analysis that is not included in this section. 

The campaign took place for four nights in September 2007, a few months before the starting 
of the turbulence monitoring runs, with both instruments, the E-ELT and lAC DIMMs, located at 



" ^http : //www ■ iac . es/proyecto/site- testing/index . php?option=com_content&task=view&id=7 
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ground level and at a distance of four meters from each other. 

Although the telescope apertures are not the same (the lAC-DIMM is 8" while the other 
is 11") the combination of telescope and CCD gave rise to a very similar pixel scales, around 
0.8 arcsec pix~' . They were empirically measured through the observation of a double star with 
known angular separation. So, similar performances of both systems were expected. The study 
was restricted to the observation of the same stars by both instruments and at the same time. The 
targets were selected to be as close as possible to the zenith. 

Both seeing time series followed a similar behavior. The differences between the seeing 
values measured by the E-ELT DIMM and those of the lAC-DIMM selt-ds ~ ^iac-dimm had a 
mean value of 0.035 arcsec, a median of 0.037 arcsec and a standard deviation of 0.099 arcsec. 

As a result, seeing values provided by the E-ELT site characterization DIMM device are 
in good agreement with those of lAC-DIMM within the measured range: from e - 0.3 arcsec 
to e - 1.1 arcsec. This is shown in Figure [H where the lAC-DIMM data are plotted versus 
those acquired by the ELT-DS DIMM. A linear fit, j = 5jc + A with the condition A = 0, 
yields a slope very close to the unity, B = 1.01 + 0.01. Unfortunately, bad seeing values (those 
worse than 1 arcsec) were scarce and therefore were not well sampled during that four nights of 
cross-calibration. 



2.1.4. Data rejection 



Raw data provided by each instru ment are validated and filte red following standard criteria. 



In the case of DIMM device, following 



Munoz-Tunon et al.iil997i. the longitudinal FWHMf and 



the transverse FWHM, seeing are compared so that only data fulfilling (Eq. |7]) are taken into 
account: 
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0.8 < ^TTTTT-^ < 1.2. 



(7) 



FWHMt 

The reason for rejection comes fror n the physics bases of the DIMM technique (see the study on 



Vernin & Muhoz-TuhonI (119951) ) and ensures the reliability of 



its uncertainties and errors made by 
the measurements. 

For what concern the MASS device, several parameters such as the removal of vignetted data 
an d the correction of M ASS overshoot are taken into account. All them are inspired by the work 



of 



Komilov et al. 



20071 The flux recorded by the MASS D channel Fo and its background signal 



Bf^ are used to check the signal-to-noise ratio. Uncertainty of the free atmosphere seeing Sf^ 
provided by the MASS software is also used as test value. Finally, the chi-square corresponding 
to the restoration of the Cjj{h) profile is also taken into account. The adopted criteria for MASS 
accepted measurements are the following: 



o Prevents too faint star or clouds: Bf^ < 3% and Fd > 100 pulses /ms. 



o Prevents too bright sky: relative Fd error < 0.03. 



o (Ts, < 0.15 arcsec. 



o Prevents bad profile restoration: from restored C^(/z) profile, < 100. 
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2.1.5. Boundary layer contribution 

The boundary layer seeing s^u here defined as the integrated turbulence between h = 5m and 
h = 500 m, is evaluated from combined MASS-DIMM observations, as follows: 

were s is the value provided by the DIMM and Sfa is gathered by integrating the MASS profiles. 
Due to noise, it may happen that doing the subtraction in (Eq. [S]) the boundary layer seeing is 
negative, and it is withdrawn from the statistics. We estimated that, so doing, any possible bias is 
almost negligible because it happened very seldom at ORM (2.8% of the whole data set), Aklim 
(3.8%) and Ventarrones (7.8%). However, at Macon this happened more often (17.6%), mainly 
during the southern winter (from August to November) coinciding with very strong wind regimes. 
It turned out that the percentage of this anomalies increases with the free atmosphere seeing, so 
rejecting those data means biasing to lower Sfa. Around the southern summertime, when the wind 
speed was lower than in winter but still higher compared to the other three sites, the percentages 
of these occurrences falls from more than 25% to around 13% of the total data acquired. See the 
next forthcoming issue of this series of papers, which will be devoted to ground meteorology, for 
a more detailed discussion. 



2.2. Complementary wind speed data 

As expressed in Sec. 11.21 and (Eq. lU, wind speed profiles are m andatory to access the 



coherence time and MASS-DIMM cannot provide these missing data. 



Travouillon et al 



2009 



wrote a long discussion about the possibility to retrieve Tq with the MASS only, according to a 



Tokovinin document but still leading to uncertainties of up to 20%. 

In order to retrieve the missing wind speed profiles necessary to solve (Eq. H]), we 
extracted them from Air Resources Laboratory (ARL) of the National Oceanic and Atmospheric 
Administration (NOAA) archivqj, from Global Data Assimilation System (GDAS) database, with 
a 1° horizontal resolution and a 3 hour temporal resolution. At ground level, we used wind speed 
given by our Automatic Weather Stations (AWS). Wind profiles have a sampling rate of 3 hours. 
In order to compare them with the MASS, DIMM and AWS databases, few assumptions were 
necessarily adopted: 1) The NOAA/ARL wind speed profiles were considered constant from 1.5 
hours before and after the V(h) time stamp. At ground level, wind speed is obtained from the AWS 
(at 2 m at Aklim and at 10 m elsewhere) to complete the vertical wind profile. 2) MASS sensitivity 
function, which approximately consists on triangular functions whose peaks are distributed on a 
2" logarithmic scale (at 0.5, 1, 2, 4, 8 and 16 km), is applied to V(h) profile in order to obtain the 
wind speed at the layers defined by MASS instrument. 3) Value of Cjj(h) corresponding to the 
first 500 m, which is not reachable by the MASS, were obtained through the combination of the 
information provided by DIMM (the total atmosphere seeing) and MASS (the free atmosphere 
seeing and the and the C^ih) from h = 500 m above). With all this information that involves 
measurements of DIMM, MASS, AWS and NOAA/ARL wind profiles, tq can be estimated (see 
(Eq. SI)). Once tq is obtained, the coherence etendue Gq can be computed (Eq. [5]). 



3. Global statistics 

In this section, the statistics of seeing (e), isoplanatic angle coherence time tq, "coherence 
etendue" Go, Fried's radius (ro), free atmosphere seeing Sfa, boundary layer seeing Sbi, at each of 



5www.ctio.edu/ atokovinin/profiler/timeconst.pdf 
'jhttp : //www . arl . noaa . gov/READYamet . php 
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the four sites: AMim, Macon, ORM and Ventarrones sites is presented. 

Data processing software has been implemented by lAC team and we made it available to 
all institutions who take care of data gathering and analysis at the different sites (see Paper I). 
Statistics of the above mentioned parameters are obtained from the whole data set and according 
to the previous remarks concerning the turbulence within the boundary layer and the wind 
profiles. From the probability distribution of each parameter the cumulative distribution and 
four percentiles: 0.05, 0.25, 0.75 and 0.95 were computed together with the mean, the standard 
deviation of the mean and the median (0.50 percentile) of the corresponding data subset. 

Many properties of these parameters might be analyzed depending on each possible purpose, 
such as the trend of the parameters along a year (Sec. H]), their typical behavior during the night 
(Sec. [5]) and their temporal stability (Sec. (6]). The global statistics only takes into account the 
valid data (those that fulfill the criteria mentioned in Sec. 12.1.41) as a whole, regardless of any 
temporal consideration. The global statistics of the parameters e, Oq, tq, Gq, ro, Sfa and Sbi are 
presented in Tables [3] and |4] over the whole observing campaign at the four sites. Both tables show 
the median, the mean, the standard deviation of the mean, cr, four percentiles, 5%, 25%, 75% and 
95%, the number of accepted data Nacc, as well as the percentage of rejected data, in accordance 
to the criteria mentioned in Sec. 12.1.41 and the total of continuous observing hours, tobs, for each 
parameter and at each site. 

Notice in Table [3] that the number of tq and Go measurements corresponding to Aklim site 
is much less than those of the other three sites. The reason for so extremely low number of data 
for these parameters in Aklim comes from the sampling of its AWS. As already mentioned in 
Sec. 12.21 the tq computation (so, also Go's) requires the wind speed at ground level acquired by 
the AWS and, while for the other three sites their sampling is a record per minute, at Aklim it is 
only one record every five minutes. This led to the aforementioned poorer statistics. 

The summary of the best achieved values of the different parameters at the different sites is 
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the following: 



During the extend of the E-ELT site campaign, the lowest median integrated seeing was 
obtained at ORM {s = 0.80 arcsec). If this particular value is compared with the results of the 
TMT site testing campaign, it is found that it is around 0. 1 arcsec hig her than those obtained 



at the TMT candidate sites showing the best total seeing statistics (see 



Skidmore et all tOO% ). 



However, in the particul ar case of the ORM, for ins tance, previous studies have proven that they 



are very similar (see e.g. 



Munoz-Tunon et al, 



(| 19971) ). The best median isoplanatic angle was very 



similar both at Ventarrones and at ORM {9q - 2 arcsec) and they are comparable to most of those 
found in the TMT candi date sites, except Mau na Kea, which benefits from a better isoplanatic 



angle (^o = 2.69 arcsec) (|Skidmore et al. 



120091) . The best (highest) median coherence time was 
measured at ORM (tq = 5.58 ms), closely followed by Ventarrones, and both of them better 
(to ^ 5 ms) than at Aklim and Macon (tq « 3.5 ms). The tq values at ORM and Ventarrones are 



also comparab 
than 5 ms (see 



e to those of the candida te sites of the TMT with highest coherence time, tq better 



Travouillon et al, 



(120091) ). Finally, the combined parameter Go defined in (Eq. [5]) 
was clearly higher at ORM (Go = 0.4 m^ ms arcsec^) and Ventarrones (Go = 0.3 m^ ms arcsec^) 
than at Macon and Aklim (Go - 0.1 m^ ms arcsec^). 

The smallest contribution of the free atmosphere measured by the MASS instrument was 
obtained again at ORM (sfa = 0.31 arcsec) although the median contribution of the boundary layer 
measured was lower at Macon (sbi = 0.51 arcsec) and at Ventarrones (ebi = 0.60 arcsec) than at 
ORM (gbi = 0.65 arcsec). The relative contribution of the ground layer and the free atmosphere 
to the total seeing at each site is shown in Table [51 In this regard, a turbulence profile showing 
a higher proportion of boundary layer turbulence, with a relatively clear free atmosphere, will 
be much more tractable for an adaptive optics syste m than one wit h , for example, strong jet 



stream-related turbulence in the tropopause (see e.g. 



Marks 



(I2OO2I) : 



Vemin & Munoz-Tunon 
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In Figures [2] to |71 the histogram as well as the cumulative distribution of seeing e, isoplanatic 
angle Oq, coherence time Tq, "coherence etendue" Go, free atmosphere seeing Sf^ and boundary 
layer seeing s^i, again, at each of the four sites are plotted. As it is well known, the conditions, 
in terms of atmospheric turbulence, are more favorable when some parameters are small: seeing 
(integrated, free atmosphere and boundary layer) or/and when they are large: isoplanatic angle, 
coherence time, Fried's radius and "coherence etendue". In the last cases (^o^ '''o, fo and Go) 
instead of estimating the cumulative distribution, the complementary cumulative distribution is 
calculated, which equals 1 minus the cumulative distribution, and is drawn in blue rather than in 
red in the plots. The four percentiles and the median are indicated by dotted lines, while the mean 
is marked with a dashed one in each of the figures. 

As a summary, the cumulative distributions of the four candidate sites were put together in a 
plot for each parameter in Figure [H From these cumulative distributions at the different sites it is 
concluded that ORM shows the best behavior in all the values e, Oq, tq and Go, closely followed 
by Ventarrones. 



4. Seasonal evolution 



Although the length of the FP6 campaign is only slightly longer than a year, the monthly 
variations of the quantities under consideration are shown in Figures |9] and \TOi in particular, 
the statistics of the parameters s, 9q, to. Go, eta and su are presented for each month along the 
whole observing campaign at the four sites. Surprisingly, the seeing seems better (lower) during 
May-Aug 2008 and Jan-Apr 2009, in both hemispheres, when one would expect an inverse trend 
depending on the hemisphere. In large databas e studies, e.g. at ORM, th e seasonal trend is more 



remarked, with summer being the best period (|Mufioz-Tufi6n et al 



1997 ). We conclude that this 



one-year is not enough for study the seasonal evolution. 
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5. Evolution during the night 

Nightly evolution of e, ^o, tq, Go, eta and Sbi at each of the four sites, over the whole campaign 
is plotted in Figures [TT] to [161 Mid-line represent the middle of the astronomical night (the middle 
point between sunset and sunrise). Due to the fact that the length of the night varies during the 
year, the beginning and the end of the figures are worse sampled than around the astronomical 
midnight (this is shown with a black curve in the plots). Every quarter of an hour, before and after 
midnight, all the data have been averaged in order to put into evidence any trend during the night. 
No clear nightly trend is visible except perhaps at Macon site where the conditions are poor at the 
beginning of the night (large seeing, low isoplanatic angle, small coherence time and thus, low 
"coherence etendue", and they get gradually better during the night until the sunset. This behavior 
is highly correlated with the evolution of wind speed at Macon, being high at the sunset and 
decreasing along the night. This issue will be discussed in more detail in the next paper dedicated 
to meteorological statistics. 

6. Temporal stability 

In Figure [IT] is drawn the "stability" of each of the following parameters, s, 9o, tq, Gq, Sf^ 
and Sbi- Stability means the average of the time interval during which a parameter, say seeing, 
remains "better" than a given value. This means, in some cases, to be lower than that value 
(integrated seeing, free atmosphere seeing and boundary layer seeing) and in other cases to be 
higher (isoplanatic angle, coherence time and "coherence etendue"). The stability plots were build 
assuming a threshold of 4 min in order to decide whether two consecutive data points of the time 
series are considered belonging to the same time interval or, on the contrary, the continuity has 
been broken. The time intervals during which each parameter remains below (or above) a given 
value were averaged. This procedure leads to the smooth curves shown in Figure [TTl 
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This concept is important to get an idea of how many time the atmospheric conditions would 
remain stable in order to carry out a particular observation that may imply a specific adaptative 
optics configuration. As an overall result, ORM site seems to exhibit higher stability than the 
three other sites, except for what concern isoplanatic angle. 

7. Summary 

The FP6 site testing campaign for the E-ELT measurements started in April 2008 and 
finished in May 2009. Four sites were characterized: Aklim (in Morocco), Macon (in Argentina), 
ORM (in Spain) and Ventarrones (in Chile). The observations were made under almost identical 
instrumentation and setup for data analysis homogeneity. The main statistical properties of the 
parameters were discussed here. The study is limited to the observations made during the FP6 
contract. 

In this sense, it is clear that a longer campaign would have been desirable in order to get rid of 
any bias produced by peculiar conditions during particular periods of time within the observations. 
A longer campaign would have made the conclusions of this study more robust. Unfortunately, 
the time spent to setup the systems in the four sites and the time constrains naturally associated 
to the ELT-DS work package resulted in a campaign slightly longer than one year. In any case, 
detailed and valuable information on these sites is provided here. Data and results from the E-ELT 
site study can be put in a more general context by making use of longer databases, when available, 
at the different sites. 

The parameters relevant for performing high angular resolution observations were obtained 
employing several instruments (MASS-DIMM and AWS installed in each of the four candidate 
sites) and the NOAA/ARL wind profile database (needed to determine the coherence time). Data 
coming from the MASS-DIMM instruments was carefully filtered by means of standard and 
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well-known criteria in order to get rid of spurious data. The DIMM instrument measurements 
were compared with a stable lAC DIMM at ORM for several night finding satisfactory correlation 
between them. 

In the case of the reference sites, ORM and Ventarrones (some kilometers away from 
Paranal Observatory), there exist large records of atmospheric turbulence conditions, although the 
discussion here is limited to the results of the FP6 contract campaign. However, it is worth noting 
that, the present work represent the first results obtained at the new sites Aklim and Macon. 

The global statistics of the high angular resolution parameters were studied as well as their 
seasonal trend, their evolution during a typical night and their time stability. 

Concerning pure statistics, the following are the ranges of the median values taken by 
the studied parameters during the campaign: the integrated seeing, e, from O.SOarcsec (ORM) 
to l.OOarcsec (Aklim); the isoplanatic angle, Oq, from 1.29arcsec (Aklim) to ~ 2arcsec 
(Ventarrones and ORM); the coherence time, tq, from 3.37 ms (Macon) to 5.58 ms (ORM); 
the coherence etendue. Go, from 0.05 m^ ms arcsec^ (Aklim) to 0.38m^msarcsec^ (ORM); the 
Fried's radius, ro, from 10.1 cm (Aklim) to 12.7 cm (ORM); the free atmosphere seeing, gfa, from 
0.31 arcsec (ORM) to 0.66 arcsec (Macon); and the boundary layer seeing £bh from 0.51 arcsec 
(Macon) to 0.77 arcsec (Aklim). Moreover, the percentages of the contribution of the boundary 
layer seeing to the total atmosphere seeing were 71% (ORM), 65% (Aklim), 50% (Ventarrones) 
and 41% (Macon). Both reference sites, ORM and Ventarrones, presented significantly higher 
median values of the global parameter Go = 0.4 m^ ms arcsec^ and Go = 0.3 m^ ms arcsec^ 
respectively, than those found at Aklim and Macon, with Go « 0.1 m^ ms arcsec^. 

The site testing campaign lasted for around a year so, although the monthly values of every 
parameter were estimated in order to study their behavior along the year, more observations would 
be obviously needed to make conclusions about seasonal trends of the sites under consideration. 
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Regarding the trend of the parameters averaged over all observation nights, it is found that 
they are very stable in ORM and Ventarrones sites. Aklim also showed good stability along the 
night, although most of the parameters seem to behave slightly better during the second half of 
the night than during the first half. A systematic variation was identified at Macon site, where 
the observing conditions are poor at the beginning and they get gradually better, being this 
phenomenon correlated with strong winds at the beginning getting weaker to sunset. 

The temporal stability of the parameters was also investigated. ORM showed generally better 
stability than the other three sites. As an example, the total seeing remained below 1 arcsec, the 
free atmosphere seeing below 0.5 arcsec, the isoplanatic angle was higher than 1.5 arcsec and 
the coherence time was higher than 5 ms for an hour, on average, at ORM (all these parameters 
considered separately; i.e. these conditions not necessarily occurring at the same time). 
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Table 1: Synthesis of DIMM data acquisition. The number of accepted data, N^cc, the percent- 
age of rejected from the total, A^rej(%), the exposure time ?exp, the median time interval between 
measurements. At, and the number of exposures per measurement, N^^p, are shown. 

DIMM 

site A^acc A^rej (%) ?exp (ms) At (s) A^exp 

Aklim 10992 21.4 5 42 200 

Macon 29723 24.4 5 100 400 

ORM 47328 11.3 5 47 200 

Ventarrones 56547 8.8 5 101 400 



Table 2: Synthesis of MASS data acquisition. The number of accepted data, A^acc, the percentage 

of rejected from the total, Nrej(%), the exposure time ?exp and the median time interval between 

measurements. At are shown. 

MASS 



site 


N 

^ 'ace 


A^rej (%) ?exp (ms) 


Ar(s) 


Aklim 


13763 


25.2 1 


63 


Macon 


94623 


5.3 1 


63 


ORM 


35962 


15.4 1 


63 


Ventarrones 


83273 


1.9 1 


63 
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Table 3:: s, Oq, tq and Go statistics obtained from April 2008 
to May 2009. The median, the mean, the standard deviation 
of the mean, four percentiles, the number of accepted data, 
the % of rejected data and the total observing time are shown. 



s (arcsec) 


med 


mean 


cr 


5% 


25% 


75% 


95% 


N 

^ 'ace 


%rej 


tob, (h) 


Aklim 


1.00 


1.09 


0.45 


0.61 


0.81 


1.28 


1.82 


10992 


21.4 


250 


Macon 


0.87 


0.91 


0.26 


0.57 


0.74 


1.05 


1.38 


29723 


24.4 


1246 


ORM 


0.80 


0.94 


0.55 


0.46 


0.62 


1.06 


2.00 


47328 


11.3 


790 


Ventarrones 


0.91 


0.96 


0.29 


0.58 


0.76 


1.10 


1.50 


56547 


8.8 


2214 


00 (arcsec) 


med 


mean 


cr 


5% 


25% 


75% 


95% 


N 

^ 'acc 


%rej 


Us (h) 


Aklim 


1.29 


1.43 


0.75 


2.66 


1.79 


0.96 


0.55 


13763 


25.2 


296 


Macon 


1.37 


1.51 


0.77 


2.81 


1.85 


1.01 


0.66 


94623 


5.3 


1705 


ORM 


1.93 


2.02 


0.77 


3.41 


2.44 


1.48 


0.87 


35962 


15.4 


669 


Ventarrones 


1.96 


2.18 


1.92 


3.97 


2.56 


1.47 


0.84 


83273 


1.9 


2626 


To (ms) 


med 


mean 


cr 


5% 


25% 


75% 


95% 


N 

^ 'ace 


%rej 


tobs (h) 


Aklim 


3.53 


5.64 


4.98 


15.87 


7.46 


2.34 


1.59 


1004 


0.0 


80 


Macon 


3.37 


3.95 


2.67 


8.58 


4.99 


2.20 


1.28 


69376 


1.5 


1227 


ORM 


5.58 


6.51 


4.14 


14.24 


8.26 


3.70 


1.82 


36802 


0.0 


619 


Ventarrones 


4.90 


5.65 


4.38 


11.77 


7.13 


3.27 


1.58 


103782 


1.0 


1807 


Go (m^ ms arcsec^) 


med 


mean 


cr 


5% 


25% 


75% 


95% 


N 

^ 'ace 


%rej 


tobs (h) 


Aklim 


0.05 


0.32 


0.66 


1.71 


0.31 


0.02 


0.01 


1004 


0.0 


80 


Macon 


0.10 


0.35 


1.20 


1.49 


0.30 


0.03 


0.01 


69376 


1.5 


1227 


ORM 


0.38 


1.02 


1.80 


4.24 


1.16 


0.11 


0.01 


36794 


0.0 


619 


Ventarrones 


0.26 


0.68 


1.42 


2.58 


0.68 


0.09 


0.01 


103782 


1.0 


1807 
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Table 4:: Statistics of ro, Sfa ^bi obtained from April 2008 
to May 2009. The median, the mean, the standard deviation 
of the mean, four percentiles, the number of accepted data, 
the % of rejected data and the total observing time are shown. 



ro (cm) 


med 


mean 


cr 


5% 


25% 


75% 


95% 


N 

^ 'ace 


%rej 


^obs (h) 


Aklim 


10.10 


10.49 


6.43 


16.49 


12.40 


7.92 


5.55 


10992 


21.4 


250 


Macon 


11.56 


11.89 


3.18 


17.74 


13.75 


9.63 


7.33 


29723 


24.4 


1246 


ORM 


12.71 


13.07 


5.04 


21.87 


16.37 


9.55 


5.06 


47328 


11.3 


790 


Ventarrones 


11.12 


11.48 


3.30 


17.41 


13.37 


9.16 


6.72 


56547 


8.8 


2214 


Sfa (arcsec) 


med 


mean 


cr 


5% 


25% 


75% 


95% 


N 

^ 'acc 


%rej 


hobs 


Aklim 


0.52 


0.63 


0.39 


0.22 


0.35 


0.77 


1.41 


13763 


25.2 


296 


Macon 


0.66 


0.79 


0.52 


0.25 


0.43 


0.98 


1.83 


94623 


5.3 


1705 


ORM 


0.31 


0.41 


0.38 


0.14 


0.22 


0.46 


0.97 


35962 


15.4 


669 


Ventarrones 


0.55 


0.65 


0.39 


0.24 


0.38 


0.79 


1.43 


83273 


1.9 


2626 


Sbi (arcsec) 


med 


mean 


cr 


5% 


25% 


75% 


95% 


N 

^ 'acc 


%rej 


hobs 


Aklim 


0.77 


0.84 


0.38 


0.34 


0.60 


1.00 


1.54 


5596 


3.8 


136 


Macon 


0.51 


0.52 


0.21 


0.20 


0.39 


0.64 


0.86 


22189 


17.6 


976 


ORM 


0.65 


0.74 


0.40 


0.31 


0.48 


0.89 


1.47 


37624 


2.8 


637 


Ventarrones 


0.60 


0.63 


0.26 


0.27 


0.46 


0.78 


1.09 


52057 


7.8 


2071 
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Table 5: Relative contribution of the boundary layer and the free atmosphere to the total seeing 
at the different sites (obtained using the median values also shown in the table). Sampling period 
from April 2008 to May 2009. 



Site 


Sbi (arcsec) 


% 


Sfa (arcsec) 


% 


total s (arcsec) 


AMim 


0.77 


65 


0.52 


34 


1.00 


Macon 


0.51 


41 


0.66 


63 


0.87 


ORM 


0.65 


71 


0.31 


21 


0.80 


Ventarrones 


0.60 


50 


0.55 


43 


0.91 
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Fig. 1 . — Seeing values measured by lAC-DIMM vs those of the DIMM part of the MASS-DIMM 
device employed for the E-ELT site characterization. The red line is a linear fit with the condition 
A=0, leading to a 1.01 slope. 




Fig. 2. — Histogram and cumulative distribution of the seeing at each of the four sites. 
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Fig. 3. — Histogram and cumulative distribution of the isoplanatic angle at each of the four sites. 
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Fig. 4. — Histogram and cumulative distribution of the coherence time at each of the four sites. 
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Fig. 5. — Histogram and cumulative distribution of the coherence etendue at each of the four sites. 




Fig. 6. — Histogram and cumulative distribution of the free atmosphere seeing at each of the four 
sites. 
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Fig. 7. — Histogram and cumulative distribution of the boundary layer seeing at each of the four 
sites. 
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e (arcsec) 6^ (arcsec) 




Fig. 8. — From top left to right bottom, cumulative distributions of the integrated seeing, the 
isoplanatic angle, the coherence time, the "coherence etendue", the free atmosphere seeing and the 
boundary layer seeing of the four candidate sites. In the cases of 6q, tq and Go the shown curves 
are the complementary cumulative distributions. Sampling period from April 2008 to May 2009. 
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Fig. 9. — Monthly statistics of e, 60, tq, and Go at the four candidate sites from April 2008 to May 
2009. Mean (diamonds), standard deviation (error bars), and median (squares) values are shown. 
Points are slightly shifted in time for clarity and months are delimited by vertical dotted lines. 
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Fig. 10. — Monthly statistics of ro, Sfa and ebi at the four sites: Aklim, Macon, ORM and Ven- 
tarrones during the whole E-ELT site characterization campaign (from April 2008 to May 2009). 
Mean (diamonds), standard deviation (error bars), and median (squares) values are shown. Points 
are slightly shifted in time for clarity and months are delimited by vertical dotted lines. 
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Fig. 1 1 . — Nightly evolution of the seeing deduced from all the nights during the whole observing 
campaign at the four sites. The median (red squares) and the mean (blue diamonds) of each time 
interval are shown together with the standard deviation of the mean (error bars) and the number of 
data (in black). 
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Fig. 12. — Nightly evolution of the isoplanatic angle deduced from all the nights during the whole 
observing campaign at the four sites. The median (red squares) and the mean (blue diamonds) of 
each time interval are shown together with the standard deviation of the mean (error bars) and the 
number of data (in black). 
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Fig. 13. — Nightly evolution of the coherence time deduced from all the nights during the whole 
observing campaign at the four sites. The median (red squares) and the mean (blue diamonds) of 
each time interval are shown together with the standard deviation of the mean (error bars) and the 
number of data (in black). 
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Fig. 14. — Nightly evolution of the coherence etendue deduced from all the nights during the whole 
observing campaign at the four sites. The median (red squares) and the mean (blue diamonds) of 
each time interval are shown together with the standard deviation of the mean (error bars) and the 
number of data (in black). 
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Fig. 15. — Nightly evolution of the free atmosphere seeing deduced from all the nights during 
the whole observing campaign at the four sites. The median (red squares) and the mean (blue 
diamonds) of each time interval are shown together with the standard deviation of the mean (error 
bars) and the number of data (in black). 



-45- 




-5 -4 -.1 -7 -1 1 ? 5 4 5 6 
U.T.-U.T. (solar midnight) (h) 






Site: Ventorrones 



Fig. 16. — Nightly evolution of the boundary layer seeing deduced from all the nights during 
the whole observing campaign at the four sites. The median (red squares) and the mean (blue 
diamonds) of each time interval are shown together with the standard deviation of the mean (error 
bars) and the number of data (in black). 
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Fig. 17. — Mean time interval during which seeings (total, free and boundary layer) are better 
(lower) than a given value, and during which isoplanatic angle, coherence time and coherence 
etendue is better (higher) than a given value. 



